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Preface 

The  Weibull  distribution  function  has  been  widely  used 
as  a statistical  tool  for  predicting  the  fracture  probabi- 
lity in  structural  materials.  Application  of  this  statis- 
tical method  to  composite  materials  addressing  nonuniform 
stress  distributions  and  the  three  parameter  form  of  the 
theory  has  been  very  limited.  I have  performed  an  evalua- 
tion of  the  three  parameter  Weibull  distribution  function 
applied  to  a composite  material  under  both  uniform  and  non- 
uniform  stress  distributions  which  includes  establishment 
of  a parameter  determination  methodology  as  well  as  the 
evaluation  of  the  theory's  ability  to  predict  fracture 
across  a variety  of  failure  modes. 

I am  deeply  indebted  to  several  individuals  who  assist- 
ed in  the  conduct  of  this  thesis.  I am  very  grateful  to  my 
thesis  advisor,  Dr.  Anthony  N.  Palazotto,  for  his  guidance, 
encouragement  and  counseling  which  aided  immensely  through- 
out this  study;  to  Dr.  James  Whitney  and  Mr.  Marvin  Knight 
of  the  Air  Force  Materials  Laboratory  for  their  support  and 
assistance  in  providing  the  experimental  data  required  to 
conduct  this  thesis.  A special  thanks  to  my  family  for 
their  patience  during  this  trying  time,  and  especially  to 
my  wife,  Jean,  for  her  many  hours  of  dedicated  typing  sup- 
port. 

Dennis  R.  Schneider 


ii 


Contents 


Page 

Preface ii 

List  of  Figures iv 

List  of  Tables viii 

List  of  Symbols ix 

Abstract xi 

I.  Introduction  1 

Background  1 

Purpose 5 

II . Theory 8 

Weibull  Distribution  Function  8 

Classical  Composite  Theory  18 

Weibull  Equation  for  Composite  Materials  . . 25 

III.  Experimental  and  Analytical  Procedures  ....  32 

Material,  Equipment  and  Test  Procedures  . . 32 

Material  Parameter  Determination  34 

IV.  Results 52 

Parameter  Determination 54 

Theoretical  Equation  Validation  66 

Alternate  Parameter  Determination  Methods  . 76 

V.  Conclusions  and  Recommendations  94 

Bibliography  98 

Appendix  Ai  Three-point  Loading  Equation  Derivation.  100 

Appendix  Bi  Experimental  Data 104 

Appendix  Ci  Data  Plots  Ill 


List  of  Figures 


Figure 


1.  Four-point  Loading  Arrangement  and  Extreme 

Fiber  Stress  Distribution  10 

2.  Stress  Distribution  - Pure  Bending 12 

3.  Third-point  Loading  Arrangement  and  Extreme 

Fiber  Stress  Distribution  .....  i4 

4.  Three-point  Loading  Stress  Distribution  ....  15 

5-  Stress  Distribution  for  Tension  17 

6.  Lamina  Axis  System  and  Cross  Section  Geometry  . 19 

7*  Typical  Stress  Distribution  - Composite  with 

Unidirectional  Ply  Orientation  ........  27 

8.  Test  Specimen  Geometry  .......  33 

9.  Parameter  Determination  Method . 38 

10.  Typical  Example  of  Weibull  Parameter  Deter- 
mination Plot 4o 

11.  Area  of  Curve  Smoothing  Sensitivity 42 

12.  Optimum  Parameter  Determination  ........  44 

13*  Comparison  of  Sample  Experimental  Data  and 
Theoretical  Cumulative  Distribution  Function 
Using  Published  Parameters  46 

14.  Comparison  of  Sample  Experimental  Data  and 
Theoretical  Cumulative  Distribution  Function 

Using  Derived  Parameters  ......  47 

15.  Error  Minimization  Process  .....  50 

16.  0-degree  Pure  Bending  Parameters 55 

17*  0-degree  Pure  Bending  Parameters  .......  $6 

18.  90-degree  Pure  Bending  Parameters  (without 

strain  gaged  specimen)  ...  57 

19*  90-degree  Pure  Bending  Parameters  (without 

strain  gaged  specimen)  . 58 


iv 


Figure  Page 

20.  90-degree  Pure  Bending  Parameters  (including 

strain  gaged  specimen)  59 

21.  90-degree  Pure  Bending  Parameters  (including 

strain  gaged  specimen)  60 

22.  Comparison  of  Experimental  Bending  Data  and 

Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen  ....  6j 

23.  Comparison  of  Experimental  Bending  Data  and 

Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen*  64 

24.  Comparison  of  Experimental  Bending  Data  and 

Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  65 

25.  Comparison  of  Experimental  Tension  Data  and 

Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen  - using  Bending  Parameters  . 69 

26.  Comparison  of  Experimental  Three -point 
Loading  Data  and  Theoretical  Cumulative 
Distribution  Function  0-degree  Specimen  - 

using  Bending  Parameters  72 

27  • Comparison  of  Experimental  Three-point 

Loading  Data  and  Theoretical  Cumulative 
Distribution  Function  90-degree  Specimen  - 
using  Bending  Parameters* 74 

28.  Comparison  of  Experimental  Three-poinx 
Loading  Data  and  Theoretical  Cumulative 
Distribution  Function  90-degree  Spe'i~en  - 

using  Bending  Parameters  75 

29.  0-degree  Tension  Parameters  78 

30.  0-degree  Tension  Parameters  79 

31.  90-degree  Tension  Parameters  80 

32.  90-degree  Tension  Parameters  81 

33.  Comparison  of  Experimental  Tension  Data  and 

Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen  82 

•including  strain  gaged  specimen  **  without  strain  gaged 

specimen 


v 


Figure  Page 

34.  Comparison  of  Experimental  Tension  Data  and 

Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen  83 

35.  Comparison  of  Experimental  Bending  Data  and 

Theoretical  Cumulative  Distribution  Function 
O-degreo  Specimen  - using  Tension  Parameters  . 86 

36.  Comparison  of  Experimental  Bending  Data  and 

Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen*  - using  Tension  Parameters  67 

37*  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  - using  Tension 
Parameters 88 

38.  Comparison  of  Experimental  Three-point 
Loading  Data  and  Theoretical  Cumulative 
Distribution  Function  ^0-degree  Specimen  - 

using  Tension  Parameters  89 

39.  Comparison  of  Experimental  Three -point 

Loading  Data  and  Theoretical  Cumulative 
Distribution  Function  O-degreo  Specimen  . . 92 

40.  Comparison  of  Experimental  Three -point 

Loading  Data  and  Theoretical  Cumulative 
Distribution  Function  00 -degree  Specimen  . . 93 

1C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
0 -degree  Specimen  112 

2C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen  113 

3C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  114 

40.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  115 

5C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  116 

•including  strain  gaged  specimen  **without  strain  gaged 

specimen 


vi 


Figure  Page 

6C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 117 

?C,  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  118 

8C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen**  119 

♦including  strain  gaged  specimen  **without  strain  gaged 

specimen 


List  of  Tables 


Summary  Correlation  Chart  53 

Weibull  Parameters  Using  Bending  5^ 

Bending  Data  Fit  Error 61 


Correlation  Results  - Bending  Parameters  . . 76 

Weibull  Parameters  using  Tension  

Tension  Data  Fit  Error  

Correlation  Results  - Tension  Parameters 

Weibull  Parameters  Using  Three-point 
Loading 

Four  Point  Loading  Data  - 0-degree  Specimen 

Four  Point  Loading  Data  - 90-degree  Specimen 

Tension  Loading  Data  - O-degree  Specimen 

Tension  Loading  Data  - 90-degree  Specimen  . 

Three  Point  Loading  Data  - O-degree  Specimen 

Three  Point  Loading  Data  - 90-degree 
Specimen  .....  . . 


En‘E22 


hk,hk-l 


Kx’VKxy 


QUt5ij 


Li3t  of  Symbols 


risk  of  rupture 

width  of  specimen  in  y-direction,  in. 

modulus  of  a unidirectional  composite  in  the 
longitudinal  and  transverse  direction  of  the 
fiber  axes,  respectively,  lb/in-. 

effective  bending  modulus  of  a laminated 
beam,  lb/in  . 

shear  modulus  of  a unidirectional  composite 
in  the  1-2  plane,  lb/in". 

thickness  of  specimen,  in. 

distance  from  the  neutral  axis  to  the  two 
outer  edges  of  the  k th  lamina  within  a 
laminated  plate,  in. 

. . 4 

beam  area  moment  of  inertia,  in  . 

lamina  number  within  a laminated  plate. 

curvature  in  a plate. 


specimen  length,  in. 

flaw  density  exponent  (shape  parameter). 

interger  less  or  equal  to  m. 

bending  moment  applied  to  a beam,  lb- in. 

bending  moment  per  unit  length  of  plate  sur- 
face, lb- in/ in. 

axial  force  per  unit  length  of  plate  surface, 
lb. 

load  on  specimen,  lb. 

lamina  in  which  a = a occurs. 

u 

terms  in  the  composite  laminate  stiffness 
matrix  in  the  material  and  load  axes, 
respectively,  lb/in41. 

dummy  variable  used  in  series  notation. 


ix 


t 


s 


probability  of  fracture. 


u,  v,  w 

V 

x,y,z 

1,2,3 


€1,€2,€12 

€ , € , € 
x y xy 


displacements  of  a specimen  in  the  x,y,z 
directions,  respectively,  in. 

volume  of  a specimen,  in^. 

load  axis  of  a specimen. 

material  property  axis  of  a specimen. 

fractional  portion  of  the  term  m. 

strains  in  the  material  axis,  in/in. 

strains  in  the  load  axis,  in/in. 

strains  in  the  mid  surface  of  a plate,  in/in. 

Poisson's  ratio  in  the  longitudinal  and 
transverse  material  directions,  respectively. 

2 

stresses  in  the  material  axis,  lb/in  . 

stresses  in  the  load  axis,  lb/in  . 

2 

stress  in  a specimen,  lb/in  . 

2 

bending  stress  in  a specimen,  lb/in  . 

mean  failure  stress  of  a distribution,  lb/in  . 

scale  parameter,  lb/in2. 

scale  parameter  for  tension  normalized  for 
volume,  lb/in  . 

scale  parameter  for  tension  normalized  for 
the  distribution,  ib/in2. 

scale  parameter  for  chree-point  loading 
normalized  for  volume,  lb.'in  . 


0 

u 


zero  probability  strength  parameter  (location 
parameter) , Ib/in  . 


* 


x 


APIT/G AE/AA/7  8M - 14 


Abstract 

XI 

This  thesis  is  an  evaluation  of  the  three  parameter 
Weibull  distribution  function  for  predicting  fracture  in  a 
composite  material  subjected  to  failure  under  both  uniform 
and  nonuniform  stress  distributions.  The  specific  forms  of 
the  three  parameter  Weibull  equations  for  these  failure 
modes  are  derived  for  a general  laminated  composite  and 
simplified  for  the  special  case  of  a unidirectional  com- 
posite. An  analysis  into  a parameter  determination  method- 
ology which  is  mathematically  reproducible  is  presented. 

The  resulting  expressions  and  methodology  are  applied  to 
experimental  bending,  tension  and  three-point  loading  fail- 
ure data  for  0-degree  and  90-degree  unidirectionally  lam- 
inated graphite -epoxy  specimer^  Weibull  parameter  sets  are 
derived  from  both  the  bending  and  tension  experimental  data. 
Each  set  is  then  used  to  evaluate  the  theory's  ability  to 
predict  the  probability  of  failure  throughout  the  three 
failure  modes  and  thereby  to  establish  a single  set  char- 
acteristic of  the  material.  v Although  a set  of  parameters 
peculiar  to  each  failure  distribution  was  obtained a char- 
acteristic set  of  values  capable  of  predicting  failure  a- 
cross  the  variety  of  failure  modes  was  not  found  for  this 
composite  material. 
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EVALUATION  OF 


THE  THREE  PARAMETER  WE I BULL  DISTRIBUTION 
FUNCTION  FOR  PREDICTING  FRACTURE  PROBABILITY 
IN  COMPOSITE  MATERIALS 


I . Introduc  tion 


Background 

The  classical  theory  of  strength  can  readily  be  3hown 
to  be  incompatible  in  material  failure  prediction  when  com- 
pared with  experimental  measurements.  This  is  most  easily 
demonstrated  by  the  dispersions  in  ultimate  strength  which 
occur  in  experimental  measurements. 

In  the  past  this  variability  of  ultimate  strength  has 
generally  been  treated  as  relatively  unimportant,  since  it 
could  easily  be  taken  into  account  through  the  U3e  of  large 
factors  of  safety.  However,  with  requirements  for  greater 
economy,  weight  restrictions,  and  increased  performance, 
significant  effort  towards  the  development  of  more  accurate 
failure  prediction  techniques  is  being  expended.  To  bridge 
the  inconsistencies  between,  the  theoretical  strength  pre- 
dictions and  experimental  data,  statistical  theories  of 
fracture  have  become  a widely  used  tool. 

The  most  widely  accepted  statistical  theory  of  frac- 
ture is  that  based  on  the  Weibull  distribution  function 
[6,2].  This  theory  is  based  on  a weakest-link  failure  con- 
cept and  relates  the  probability  of  fracture  to  the  actual 
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stress  observed  at  fracture  [2].  To  fit  the  physical  facts 
ot  the  experimental  variability,  the  theory  has  been  cor- 
rected to  a scmi-empirical  distribution  based  on  experi- 
mental data  [2].  Two  basic  criteria  of  fracture;  size  and 
normal  tensile  stress,  are  used  in  the  theory,  and  it  is 
postulated  that  failure  in  an  isotropic,  homogeneous  mate- 
rial is  fully  described  by  three,  material  dependent  param- 
eters; the  zero  probability  strength  (location  parameter), 
the  flaw  density  exponent  (shape  parameter)  and  a scale 
parameter  [?].  Within  the  validity  of  these  criteria,  the 
theory  can  describe  failure  for  any  type  of  stress  distri- 
bution [6]. 

Use  of  this  theory  has  been  applied  to  a large  variety 
of  materials,  from  ceramics  to  laminated  composites.  Daniel 
and  Weii  [6]  investigated  the  effects  of  a stress  gradient 
on  the  fracture  characteristics  of  brittle  materials;  Bortz 
and  Weii  [2]  investigated  the  applicability  of  Weibull *s 
semi -empirical  distribution  function  to  ceramic  oxides, 
deriving  the  Weibull  material  parameters  under  various  con- 
dition of  test  temperature,  heat  treatment,  surface  condi- 
tions as  well  a3  specimen  size;  Kaminski  [ii]  studied  the 
effects  of  specimen  geometry  on  the  predictability  of  the 
strength  of  a composite  material;  and  Knight  and  Hahn  [i2] 
performed  experiments  on  randomly -distributed  short  fiber 
composites,  looking  at  the  statistical  characteristics  of 
the  experimental  data  by  U3C  of  the  Weibull  distribution 
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function. 


Most  of  the  efforts  investigating  fracture  using  stat- 
istics have  dealt  with  homogeneous  isotropic  materials,  with 
a lesser  empnasis  towards  composite  materials.  However, 
with  present  day  interest,  especially  within  the  aerospace 
industry,  being  placed  on  such  factors  as  weight,  strength, 
and  performance,  the  composite  material  has  taken  on  an 
important  role  in  the  development  of  lighter-weight  high- 
strength  aerospace  structures.  Even  though  modern  manufac- 
turing procedures  of  this  material  allow  for  carefully  con- 
trolled manufacturing  conditions,  its  precise  breaking 
strength  cannot  be  accurately  predicted  L9 3*  The  variabil- 
ity of  breaking  strengths  in  a composite  is  identical  in 
nature  to  that  found  for  other  materials,  even  when  speci- 
mer& are  tested  under  assumed  identical  test  conditions.  It 
is  this  variability  among  a controlled  population  which  nec- 
essitates accurate  statistical  approaches  to  the  fracture 
process. 

The  Weibull  distribution  function  can  be  applied  to 
experimental  data  in  either  a two  parameter  or  three  param- 
eter form.  Each  characterization  has  been  investigated 
using  both  uniform  and  nonuniform  stress  distributions. 

A discussion  of  the  parametric  equations  and  related 
functions  is  provided  in  Section  II  of  this  thesis,  however, 
it  is  important  at  this  point  to  mention  some  previous 
efforts  relative  to  the  Weibull  theory. 
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Daniel  and  Weil  [6]  addressed  in  detail  the  three 
parameter  form  of  the  Weibull  equation  for  both  uniform  and 
nonuniform  stress  distributions.  Columbia  Resin  (CR-39) 
was  used  as  the  experimental  material.  The  theoretical 
derivations  for  different  stress  fields  are  presented  based 
on  elementary  beam  theory  applied  to  a prismatic  beam  under 
four-point  loading.  The  resulting  set  of  statistical  equa- 
tions provides  a means  with  which  the  probability  of  frac- 
ture can  be  predicted  for  a variety  of  loading  conditions 
and  8 tress  distributions  assuming  the  three  material  param- 
eters are  known.  The  hypothesis  in  their  development  is 
that  the  three  parameters  are  in  fact  "material”  constants. 
This  is  substantiated  by  graphically  confirming  a close 
correlation  between  theoretical  and  experimental  data.  Al- 
though a set  of  three  material  parameters  was  established 
and  verified  for  one  set  of  experimental  data,  no  verifica- 
tion of  this  set  of  parameters  with  data  from  different 
stress  distributions  was  provided.  Similar  efforts  have 
since  been  completed  using  composite  materials,  however, 
detailed  theoretical  validation  is  very  limited. 

Daniel  and  Weil  also  developed  analytical  means  for 
determining  the  unknown  material  parameters  from  a given 
set  of  experimental  data.  Their  method  involved  the  simul- 
taneous solution  of  a system  of  three  equations  with  three 
unknowns.  However,  due  to  computational  difficulties  which 
were  experienced,  the  graphic  method  devised  by  Weibull  [16] 
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proved  to  be  much  easier  for  parameter  determination. 

Since  the  composite  material  provides  a material  with 
unique  characteristics  when  compared  to  conventional  struc- 
tural materials,  the  applicability  of  the  approach  develop- 
ed by  Daniel  and  Weil  is  of  significant  interest.  Further- 
more, the  verification  of  whether  a given  set  of  material 
parameters  will  allow  the  prediction  of  fracture  probability 
throughout  a variety  of  loading  conditions  would  confirm 
that  the  parameters  are  in  fact  material  related.  More 
important,  this  evidence  could  provide  a significant  cost 
savings  in  terms  of  the  number  tests  required  to  determine 
general  fracture  characteristics  of  a material. 

Purpose 

There  are  four  separate  purposes  to  be  investigated  in 
this  thesis. 

The  first  purpose  is  to  establish  the  basic  forms  of 
the  Weibull  equation  for  a composite  lamina  under  uniform 
and  nonuniform  stress  distributions  using  the  apprr  ch  de- 
veloped by  Daniel  and  Weil  for  a homogeneous  isotropic  beam 
[6].  Classical  laminate  theory  is  expanded  in  the  Weibull 
equation  form  for  bending,  tension  and  three-point  loading 
conditions.  The  equations  are  derived  using  general  com- 
posite material  notation  and  then  reduced  to  a simplified, 
single  lamina  form,  for  the  specific  specimen  being  used  in 
this  study.  Composites  with  only  unidirectional  ply  orien- 
tations (0 -degree  and  90-degree)  are  investigated. 
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The  second  purpose  is  to  obtain  experimental  fracture 
data  for  graphite -epoxy  to  allow  validation  of  the  above 
equations.  The  material  parameters  are  derived  from  experi- 
mental pure  bending  data  gathered  to  form  the  characteristic 
fracture  distribution  of  the  material.  Twenty-five  tests 
will  form  this  distribution.  Both  the  analytical  methods 
described  in  [7]  as  well  as  the  trial  and  error  graphical 
method  recommended  by  Weibull  [16]  are  used  to  determine 
the  material  parameters.  Procedures  to  permit  reproducibil- 
ity of  the  parameters  and  parameter  associated  sensitivities 
are  addressed. 

The  third  purpose  of  this  thesis  is  to  investigate  the 
extendability  of  the  three  parameters  found  from  bending 
tests  to  other  failure  modes.  This  will  investigate  the 
hypothesis  that  the  parameters  are  in  fact  material  related 
numbers.  A minimum  of  2 5 test3  for  both  tension  and  three- 
point  loading  will  establish  the  fracture  distributions  for 
each.  The  theoretical  equations  developed  for  the  first 
purpose  of  this  study  will  be  solved  using  the  three  param- 
eters found  from  the  bending  test  distribution.  The  crite- 
ria for  conformance  is  the  closeness  of  fit  between  each 
experimental  distribution  and  theoretical  prediction.  Meth- 
ods of  obtaining  additional  accuracy  of  fit  will  also  be 
evaluated . 

The  fourth  purpose  is  to  compare  the  three  parameter 
and  two  parameter  Weibull  function  in  prediction  accuracy. 
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The  two  parameter  model  also  considers  nonuniform  stress, 
however  the  third  parameter,  the  location  parameter,  is 
assumed  to  be  zero.  Factors  in  comparison  include  the  math- 
ematical difficulty  of  establishing  the  three  parameter 
model  and  variations  in  fit  accuracy  when  compared  to  the 
two  parameter  model.  In  addition,  suitability  of  each  model 
across  the  variety  of  failure  modes  is  to  be  addressed. 
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II . Theory 


The  theory  development  for  this  investigation  is  pre- 
sented in  three  sections.  The  first  two  sections  will  pro- 
vide a brief  background  into  the  basic  Weibull  theory  and 
classical  composite  theory.  The  third  section  will  provide 
the  derivations  and  equations  used  in  this  thesis,  which 
are  a combination  of  the  theories  described  in  the  first 
two  sections. 

Weibull  Distribution  Function 

This  section  is  intended  to  provide  a background  of 
the  theoretical  considerations  of  the  Weibull  theory  for 
the  prediction  of  fracture  [6,16]  as  well  as  its  extensi'. 
to  a variety  of  stress  distributions. 

The  Weibull  theory  [lo]  is  based  on  the  concept  tha'u 
brittle  materials  contain  a large  number  of  flaws,  which 
tend  to  lower  the  fracture  stress  of  a material  below  the 
theoretical  rupture  stress  [2].  These  flaws  are  assumed 
to  be  of  a random  size  and  distribution  throughout  the  body 
and  to  be  the  cause  of  the  scatter  which  is  found  in  experi 
mental  fracture  data.  Two  basic  critera  of  failure  used  in 
the  Weibull  distribution  function  are  size  and  normal  ten- 
sile stress.  For  a uniaxial  stress  field  in  a homogeneous 
isotropic  material,  the  probability  of  fracture  for  a vol- 
umetric flaw  distribution  is  given  by 


where 


B - j^fdv  121 

and  B is  the  risk  of  rupture 

au  is  the  zero  probability  strength  (location 
parameter) 

Oq  is  the  scale  parameter 

m is  the  flaw  density  exponent  (shape  param- 
eter) 

Equation  (1)  shows  that  no  special  allowance  is  made 
for  nonuniformity  of  stress.  Each  elemental  volume  with  a 
normal  tensile  stress  is  added  to  form  the  overall  risk  of 
rupture  for  the  specimen.  To  address  the  question  of  non- 
uniformity  [6],  one  can  look  at  specific  case  derivations 
using  a prismatic  beam  under  desired  loading  conditions  or 
stress  distributions.  The  three  loading  conditions  which 
are  investigated  in  this  thesis  are*  pure  bending,  three- 
point  loading  and  tension. 

The  development  of  equations  accounting  for  tensile 
stress  gradients  can  be  accomplished  using  elementary  beam 
theory.  Figure  1 shows  a prismatic  beam  under  four-point 
loading  and  the  resulting  extreme  fiber  stress  distribution. 
The  section  between  the  minor  span  is  assumed  to  be  in  pure 
bending,  while  the  section  between  the  minor  and  major  span 
has  the  type  of  stress  distribution  found  in  a three-point 
loading  condition.  These  sections  can  now  be  addressed 
separately. 
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For  a pure  bending  failure,  let  (k)  in  Figure  1 be 
equal  to  infinity,  thereby  placing  the  total  length  of  the 
beam  under  pure  bending.  Using  elementary  beam  theory,  the 
stress  distribution  within  the  beam  is  described  by 

<r  * z o) 

this  assumption  can  be  realized  in  a four-point  bending 
test  using  only  failures  occurring  between  the  minor  span  as 
the  data  samples. 

The  risk  of  rupture.  Equation  (2)  can  now  be  rewritten 
in  integral  form  as  follows 


where 


Zu"  zcr; 


(4) 


(5) 


Equation  (5)  is  necessary  since  the  three  parameter  form  of 
the  Weibull  equation  is  being  used,  thereby  making  the  as- 
sumption that  a stress  value  exists  below  which  the  proba- 
bility of  fracture  is  zero.  A physical  interpretation  of 
Zy  is  shown  in  Figure  2 in  which  the  shaded  area  represents 
the  region  where  the  zero  probability  failure  parameter, 
(ou)»  is  exceeded.  The  stress  gradient  for  pure  bending  is 
as  shown  for  every  point  along  the  x - axis.  Integrating 
Equation  (4)  results  in  the  risk  of  rupture  equation  for 
pure  bending 


i 

i 
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Pig,  2.  Stress  Distribution  - Pure  Bending 
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(6) 


B- 


V 


Z(m+i)  Cfc  c" 

where  V is  a dimensionless  number  expressing  the  quantity 
of  unit  volumes  subjected  to  a uniform  tensile  stress  £2]. 
Using  Equation  (6)  in  Equation  (l)  the  probability  of  frac- 
ture for  a beam  subjected  to  pure  bending  can  be  found. 

The  derivation  of  the  risk  of  rupture  for  a center  or 
three-point  loading  condition  requires  that  (k)  in  Figure  1 
be  equal  to  two.  This  is  representative  of  a beam  with  a 
varying  moment  as  shown  in  Figure  3.  According  to  elemen- 
tary beam  theory,  the  stress  distribution  within  the  beam 
is  described  by 


(T=  ^ *Z 


0 i y i f 


(?) 


Equation  (2)  can  now  be  rewritten  using  this  definition  and 
the  following  integral  limits 


B=^ 


where 


4 / V So 


v/  Cu  H L 

Xu"TKT 


dv 


(8) 


(9) 


and  2^  is  given  in  Equation  ( 5) 


The  limit  expressed  by 
Equation  (9)  represents  the  point  in  the  x - direction  at 


which  the  stress  reaches  the  value  3^.  A physical  inter- 
pretation of  xu  is  shown  in  Figure  4.  The  superimposed 
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Beam  Surface 
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Three-point  Loading  Stress  Distribution 


stress  curves  show  how  the  stress  gradient  and  extreme 
fiber  stress  varies  with  respect  to  the  x and  y axis  for 
this  loading  condition.  Using  simple  proportions  it  can  be 
shown  that 


ft/Z  = cVl»i  (10) 

from  which  Equation  (9)  is  derived.  Integrating  Equation 
(8)  results  in  the  risk  of  rupture  equation  for  a third- 
point  loading  condition 


m ■ m (I?) 

and  (m)  is  the  largest  integer  less  than  or  equal  to  (m)  and 
{*)  is  the  decimal  fractions  of  (m).  The  term  ( r)  is  intro- 
duced as  a dummy  variable  for  the  summation  term. 

Due  to  the  complexity  of  the  mathematics  involved  in 
the  derivation  of  Equation  (11),  a detailed  derivation  is 
provided  in  Appendix  A. 

The  derivation  for  the  risk  of  rupture  under  tension 
loading  is  the  simplest  of  the  three  loading  conditions 
being  investigated.  Figure  5 3hows  that  the  stress  distri- 
bution for  this  loading  results  in  a uniform  stress  over  the 
entire  thickness  of  the  beam  and  therefore  is  expressed  as 
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C * Gl 


(13) 


The  risk  of  rupture  equation  therefore  becomes 


(14) 


Integrating  Equation  (14)  results  in  the  risk  of  rupture 
for  tension  loading 


B - v 


(15) 


Once  the  failure  distribution  of  a material  has  been 
established,  another  characteristic  of  interest  is  th^  dis- 
tribution mean.  Once  the  three  material  parameters  m,  au 
and  o0  have  been  established,  the  mean  failure  stress  cm 
can  be  determined  and  used  as  a further  check  on  the  use  of 
the  Weibull  theory  for  a given  set  of  data  [2], 

oo 

Cm  ” Ca  + f 6 

'(m 

Equations  (6),  (11)  and  (15)  represent  the  basic  ap- 
proach and  format  of  the  risk  of  rupture  expressions  which 
will  be  extended  to  composite  materials. 


Classica-i.  Composite  Theory 

This  section  outlines  some  of  the  basic  laminated  plate 
relationships  (for  background  see  [l,3,10].  Figure  6 
indicates  the  axis  system  relationships  which  are  used  in 
the  derivations,  as  well  as  some  of  the  geometry  of  a 
laminated  composite. 
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Fig.  6.  Lamina  Axis  System  and  Cross  Section  Geometry 
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The  stresses  in  an  orthotopic  composite  plate,  or  lam- 


ina, can  be  calculated  with  the  following  relationships 


Qu 

Qj. 

CU 

f 0 

Qu 

Qa 

Q* 

4 

! f4 

Q(4 

Oh 

Qu 

K 

/•] 

(17) 


where  (Q.  .)  is  referred  to  as  the  lamina  stiffness  matrix 

^ J 

and  whose  components  are 


Qu  ^ Em  / Cl  - V.iVi.) 

Qzi  = Eu /(  I ” Va  V*<) 

Q44  = 6*2.  ( 13) 

Qi2  = VaEn  / C l ' VizVu)  c Vn  Eh  /Cl  “ V12V21) 

Q«*  * Q24  * 0 


and  (k)  the  k th  lamina  within  a laminated  composite.  The 
matrix  of  Equation  (1?)  provides  the  equations  for  finding 
the  stress  due  to  loading  the  material  axis.  If  the  stress 


is  desired  due  to  loading  in  other  than  the  material  axis, 
(x,y,z  in  Figure  6),  rotated  some  angle  0 off  the  material 


axis,  the  stress  is  found  using  the  following  relationship 


Q.i 

Qu 


Qu  Q it 
Qa.  Qi* 
Q*  0 


64 


< 

K 

(19) 


where 


Q.,  = Qu  Gs4e  - 2 (QtfZQd  Si^ecos'e  ♦ Gu  5"4© 

On-  0.  v© ♦ 2 C0a+2o<t)  $ m*e Cos  a + c» 

0t6r  S«20Cos2(9  + (rs40) 

Q.i-(G,.4Q2i-4Qtt)  S.*leG*2©* 

5%w  (Qm“  Qa.’ZQu)  Sm0Cos30  •+  (Gii“Q2i^206t)  5»h50Cos© 
5*6* CQu“Qrt."2G6t)  Sw30Cos0  + CQ^Cta^Cta)  Sm  0Cos30 


The  convention  used  for  positive  rotation  is  shown  in  Figure 

6. 
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The  strain  components  can  be  obtained  from  the  follow- 


ing relationships 

r 
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' o “< 
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K? 

K>y 


(21) 


where  the  (°)  superscript  refers  to  the  midplane  strains  of 
a laminated  composite.  The  components  of  Equation  (21)  are 
as  follows 


(22) 


and 


•\ 

l 

J i 

5 i 

^ y •»£ 

J‘V 


(23) 


J c2.^/Wy 

where  u,  v,  and  w are  displacements  in  the  x,  y,  and  z 
direction  respectively*  If  one  uses  the  relationships  of 
Equations  (21),  (22)  and  (23),  Equation  (19)  can  be  written 
as 


WrMkM  12*> 

For  convenience,  a simpler  system  of  forces  and  mo- 
ments is  now  introduced  which  are  defined  as  force  resultants 
and  moment  resultants.  The  force  resultant  (N  ) is  defined 
as  a force  per  unit  width  and  the  moment  resultant  (M  ) as 
a moment  per  unit  length.  In  general  then,  these  values 
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can  be  given  in  terms  of  the  stress  vector  as  follows 


where 

i = ( v , v , yy) 


By  rewriting  these  equations  * a shorter  notation  can  be 


obtained 


where 


and 


♦[>]{*} 

i-v.) 

^g*^L  LqUk-k.) 


(27) 


(28) 


[A]  - extensional  stiffness 

[B]  - coupling  stiffness 

[D]  - bending  stiffness 

- h^ _^)  represent  the  individual  lamina  thicknesses. 

When  the  laminate  is  symmetric  in  bo th  geometry  and 
material  properties  about  the  raid-plane,  Equation  (27) sim- 
plifies considerably.  Due  to  the  symmetry,  the  thickness 
terms  cause  the  coupling  stiffness  (B. .)  to  go  to  zero  [10]. 

Por  the  uncoupled  axial  and  moment  loads,  Equation  (27) 
w.-^n  be  written  as 


and  (^ 
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(29) 


M-MM 
M’MM 

Whitney,  Browning  and  Mair  [17 J further  simplified  the 
moment  resultant  equation  as  follows*  since  from  Equation 
(29) 


M •[»]>} 


(30) 


Thus,  for  a pure  bending  moment  this  can  be  written  as 

Mx 

0 

o 

-J  % 

M is  the  moment  resultant  defined  by 
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(32) 


If  one  assumes  that  the  specimen  has  a high  length- to-width 
ratio  then 

w * W(x) 

and 


dlW/ 


where 


M 

n 


ax* 

On 


(33) 


which  is  of  the  same  form  as  classical  beam  theory  with  the 
isotropic  modulus  (E)  replaced  by  an  effective  bending  mod- 
ulus (E) . Stress  can  now  be  obtained  as 


M„-[a]„z{K} 


(34) 
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or 


C nek  * Z fiK^J 


n 

I 
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where 


(Tyk  "I  f2ky 

CTik  “ Z f 3h y 

{fi}k*  [a.] 

i'  1,2,5 


(35) 


(36) 


►v»«  ‘ ,l,b 

Q<k\*  0 WA. 

S»*  - Dim  In?/ K 2. 

Equation  (36)  then  represents  an  equation  applicable  to  a 
general  class  of  symmetric  laminates,  derived  by  consider- 
ing a beam  as  a special  case  of  a laminated  plate. 

The  axial  stress  equation,  when  considering  only  the 
normal  force  resultant  N^,  can  be  derived  tt&  follows 


10  * M" 


where 


N*  • | 


(37) 


(38) 


Since  symmetric  specimens  will  be  used  the  equation  for 
tension  only  is 

{0*{£}rWNx  (39) 

the  normal  stress  can  then  be  expressed  as 


(Txic  * Cfi*K  £ * * Q\l k €y  + QtfvCxy 

S ^ Q Art.  * 0 A~*j  N* 


(40) 
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Weibull  Equation  for  Composite  Materials 


This  section  deals  with  the  extension  of  the  Weibull 
distribution  function  forms,  presented  for  a beam  in  Equa- 
tions (6),  'll)  and  (15).  to  composite  materials.  Deriva- 
tions are  based  on  the  composite  theory  equation  forms  pre- 
sented in  the  previous  section.  The  final  equations  are 
then  presented  both  in  composite  notation  as  well  as  in  a 
simplified  form  applicable  to  the  specific  experimental 
specimens  being  analyzed  in  this  study. 

Since  the  Weibull  theory  is  based  on  fracture  duo  to 
normal  tensile  stress  in  materials  containing  random  flaw 
distributions,  the  composite  is  also  assumed  to  contain 
volumetric  flaws  and  to  be  subject  to  normal  tensile  stress 
fracture.  The  derivation  of  the  composite  bending  equations 
can  be  modified  with  this  assumption,  and  Equation  (35)  re- 
duces to 


or 

u **  s 2.  D»»  * Da  4 D*J  (42) 

if  we  now  define 

C*  * [Q..,D’.  *Q-:k  da  ♦ 5..kD1]  ^ t*3> 


for  simplified  notation,  then  the  stress  in  a laminate  sub- 
jected to  bending  can  be  expressed  as 


<u  * Z Ck 


(44) 
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where  (2)  is  the  distance  from  the  neutral  axis. 

The  laminate  specimen  used  in  this  investigation  will 
consist  of  only  unidirectional  ply  orientations  and  there- 
fore act  like  a single  lamina  with  respect  to  stress  dis- 
tribution. Equations  (45)  through  (49)  provide  some  basic 
definitions  for  later  theory  development  and  are  graphically 
presented  for  a unidirectional  laminate  in  Figure  7 


(45) 


where  the  subscript  (p)  refers  to  the  p-th  lamina  which  is 

the  lamina  in  which  0 = a . The  stress  in  the  p-th  lamina 

u r 

is  defined  as 


(T  * CVZ* 

[Gv  D»* 4 QiipDu  ♦ GU?  D*] 

M 

b 

(46) 

also 

(Juf  9 Cp  h<rH 

(4?) 

Gh  b C* 

(48) 

(Tp  9 Cp  hp 

(49) 

With  the  use  of  these  equations  the  evaluation  of  the 

ri3k 

of  rupture  becomes 

8- 

0 0 w 1 

0 0 V-i 

1 w , 

I dv 

(50) 

where  the  first  triple  integral  represents  the  number  of 
unit  volumes  in  the  lamina  in  which  the  stress  partially 
exceeds  ou«  The  second  triple  integral  represents  the  sum- 
mation of  this  quantity  for  all  other  lamina  in  which  the 
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with  Unidirectional  Fly  Orientation 


stress  exceeds  a^.  Subscripts  are  provided  to  au  to  ac- 
count for  variations  of  this  value  in  multi-directionally 
oriented  specimens  and  as  such  make  this  term  ply  dependent 
The  equation  in  this  form  represents  the  risk  of  rupture 
for  pure  bending  of  a general  laminated  composite.  For  l 
unidirectionally  laminated  composite,  as  that  being  used  in 
this  thesis,  au  is  assumed  to  be  identical  for  all  plies 
and  consequently  subscripts  are  dropped  from  subsequent 
equations.  Integrating  this  equation  results  in  the  expres 
sion  for  the  risk  of  rupture  for  a unidirectionally  laminat 
ed  composite  subjected  to  pure  bending. 


.Mui 


f rrckh.-(Tui">ti-rah^-6~u)"'<: 
.1  Cp  J Cp 


(5D 


This  equation  can  be  simplified  further  by  use  of  the  char- 
acteristics peculiar  to  the  specimens  being  incorporated 
into  this  thesis.  For  unidirectional  fiber  orientations 


Cp 


C 


(52) 


because  the  terms  in  the  lamina  stiffness  matrix  Q. . are 

• 0 

identical  for  each  lamina.  This  reduces  Equation  (51)  to 
rs  Lb  (chn-<m)  f Chn~  (53) 

B ■ JmVf  C \~G~I 

Equation  (53)  can  also  be  written  in  terms  of  the  extreme 
fiber  stress  as  was  Equation  (6).  The  stress  distributions 
for  these  composite  specimens  are  linear,  and  therefore  can 
be  expressed  as 
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since 

r cfe)=  ^ 

therefore 

n Lb  /fo-(TuW<r«-(r^w 
D=  ISiTi)  (— ) i-0T7 

where 

C=  (Tn(2/h) 

now 
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^ 2(m+ 1)  y Cn  / y Co  j 

where 

Lbb  =V 


(55) 


(56) 

(57) 


(58) 


Equation  (58)  is  identic rl  to  Equation  (6)  as  is  expected 
for  the  type  of  composite  being  considered,  and  therefore 
represents  the  risk  of  rupture  equation  for  pure  bending 
loading  tests  in  the  thesis.  A similar  approach  is  used  in 
the  derivation  of  the  risk  of  rupture  equation  for  three- 
point  loading.  Some  of  the  previously  defined  terms  must, 
however,  be  modified  to  allow  for  a varying  moment.  The 
moment  equation  for  three-point  loading  is 

m - fyy  (59) 

where  (P)  is  the  applied  load.  Equation  (42)  can  now  be 
written 
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(60) 


defining 


then 


(Tk  * GZX 


(61) 


(62) 


since  the  specimen^ are  unidirectional,  the  lamina  subscript 
can  be  dropped.  Figure  3b  shows  the  extreme  fiber  stress 
distribution,  and  from  this  we  define 


Gb  (beam)*  (Jn (composite  ) 


r*  _bnL 

C 2 


(63) 


or 

c'.  (64) 

hn  L 

where  o is  the  maximum  extreme  fiber  stress,  and  h,  is  de- 
n n 

fined  in  Figure  ?.  Using  Equation  (62)  in  Equation  (2)  re- 
sults in  the  expression 


which  when  integrated  takes  the  form 


(66) 


which  can  be  put  in  terms  of  the  maximum  extreme  fiber 
stress  an  by  use  of  Equation  (64).  Using  this  definition 
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the  final  risk  of  rupture  for  a composite  under  three-point 
loading  then  becomes  identical  to  Equation  (11). 

The  risk  of  rupture  for  a composite  under  tension  is 
derived  from  Equations  (25)  and  (27).  In  this  case  only 
the  tension  loading  and  extensional  stiffness  is  considered. 
Using  Equations  (2)  and  (40)  the  integral  equation  for  the 
risk  of  rupture  under  tension  is 


which  integrates  to 

B = Lb£  ( (68) 

K:' 

If  one  incorporates  the  unidirectional  properties  of  the 


specimen  as  related  to  this  thesis  in  addition  to 
(J"*k  : ft  '• 


(69) 


Equation  (68)  reduces  to  Equation  (15)* 

This  study  will  utilize  Equations  (15)»  (58)  and  (66) 
in  predicting  the  risk  of  rupture  for  graphite-epoxy  under 
the  loading  conditions  for  which  each  of  the  equations  was 
derived,  and  then  use  these  values  in  establishing  the  prob- 
ability of  fracture  from  Equation  (1). 
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Ill . Experimental  and  Analytical  Procedures 


This  section  provides  an  overview  of  the  methodology 
developed  to  obtain  the  material  parameters  in  the  Weibull 
distribution  function.  A description  of  experimental  equip- 
ment and  procedures  as  well  as  the  analytical  process  used 
in  manipulating  the  experimental  data  is  provided.  Also  pos- 
tulated in  this  section  is  some  analysis  theory  directed  at 
data  evaluation  and  interpretation.  A sensitivity  analysis 
is  provided  in  an  attempt  to  show  the  dependency  of  the  re- 
sulting theoretical  equations  on  the  specific  experimental 
data  being  analyzed. 

Material,  Equipment  and  Test  Procedures 

Graphite -epoxy  specimens  are  used  exclusively  in  the 
experiments.  Specimens  were  constructed  of  both  sixteen 
and  eight  ply  graphite  fibers  impregnated  in  epoxy.  Fiber 
orientations  were  unidirectional  within  each  specimen.  The 
specimens  were  cut  from  a single  24  inch  by  24  inch  sheet, 
thereby  obtaining  both  0-degree  and  90-degree  oriented  test 
samples  from  the  same  layup.  A description  of  the  specimens 
is  provided  in  Figure  8.  Two  laminate  thicknesses  are 
used  due  to  limitations  on  specimen  sizei  sixteen  ply  spec- 
imen for  both  four-point  and  three-point  loading  tests  and 
eight  ply  specimens  for  the  tension  tests.  The  lower  ply 
specimens  were  chosen  for  tension  so  the  necessary  fracture 
loads  would  not  deform  or  separate  the  end- tabs  from  the 
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.08" 


Typical  16  Ply  90  & 0 degree  Specimen 
(Pure  Bending  & 3-point  Loading) 


. 04” 


Typical  8 Ply  90-degree  Tension  Specimen 


04m 


Typical  8 Ply  0-degree  Tension  Specimen 


Fig.  8.  Test  Specimen  Geometry 
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specimen,  thereby  inducing  measurement  inaccuracies. 

A standard  Instron  testing  machine  was  used  to  conduct 
all  testing.  For  each  four- point  and  three-point  loading 
test,  the  load  was  applied  at  .05  in/min.  The  tension  tests 
were  conducted  with  a stroke  rate  of  .02  in/min.  The  appli- 
cation of  loading  continued  until  castastrophic  failure  or 
load  resistance  went  to  zero.  A chart  recorder  recorded 
deflection  versus  load  for  each  specimen.  Approximately  20$ 
of  the  specimens  were  equipped  with  strain  gages.  For  these 
specimens  strain  versus  load  was  recorded. 

Accurate  volumes  were  obtained  for  each  individual 
specimen  by  conducting,  individually,  three  measurements  of 
thickness  and  width.  These  measurements  were  then  averaged 
to  obtain  an  average  width  and  thickness  per  specimen.  All 
testing  followed  in  accordance  with  standard  test  proce- 
dures. In  order  to  obtain  a region  of  pure  bending  within 
the  four-point  loading  tests,  the  cross-sectional  dimensions 
were  controlled  to  a span  to  depth  ratio  of  32  to  1 for  each 
specimen.  A minimum  of  25  specimens  were  tested  for  each 
ply  orientation  and  loading  condition. 

The  succeeding  paragraphs  provide  the  steps  followed  by 
the  author  in  arriving  at  the  final  procedures  of  deriving 
the  material  parameters  using  published  information.  With 
the  accumulated  experience  obtained  in  the  methodology  pro- 
cess, a final  technique  is  developed  and  presented. 

Material  Parameter  Determination 

Once  the  basic  equations  for  the  Weibull  distribution 


function  are  established,  experimental  data  is  gathered  to 


validate  the  theory.  The  problem  that  now  arises  is  the 


development  of  those  values  of  the  material  parameters  which 


make  the  theoretical  probability  of  fracture  curve  fit  the 


experimental  data  best. 


Daniel  and  Weil  [6]  used  as  a suitable  criterion  the 


analytical  method  of  minimization  of  the  sum  of  the  mean 


squares  differences.  In  their  procedure,  experimental  data 


obtained  from  pure  bending  tests  was  used  in  the  parameter 


determination. 


The  procedure  combines  Equation  (1)  and  (6).  Taking 
the  double  log  of  Equation  ( 1)  yields 

Y"  --  LnU  yr^  = Ln  B : In  yr  - LnO*1)  (?0) 
+ (lY Hi)  Ln  - In  Gn  * 


The  corresponding  (estimated  value  for  the  probability  of 


fracture  obtained  experimentally  is 

< -JL- 

•in-  M+  I 

from  which 


(71) 


(72) 


where  N = total  number  of  specimen/  tested 

n = specimen  number  when  specimen  are  arranged  in 


ascending  order  of  stress 


The  least  squares  method  requires 

w 


£ (Yn ■ y*)1--  MINIMUM 


(73) 


The  necessary  conditions  for  thu  minimization  also  require 


that  the  derivatives  of  Equation  (73)  when  taken  with 


respect  to  ou»  0Q  and  m to  be  set  equal  to  zero. 
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This  then  results  in  three  equations  with  three  unknowns 
which  can  be  solved  simultaneously.  Convergence  to  zero 
results  in  the  determination  of  a unique  parameter  set. 

Although  mathematically  correct,  convergence  of  these 
equations  was  found  difficult  to  obtain.  A possible  reason 
being  their  highly  non-linear  nature  [2]. 

Another  more  widely  used  approach  is  one  suggested  by 
Weibull  [6].  This  method  is  a trial  and  error  procedure  in 
which  the  experimental  data  is  linearized.  The  recommended 
procedure  is  to  rearrange  equation  (70)  into  the  form  of 
the  equation  of  a straight  line 

y = na  x +b  (74) 

this  is  accomplished  as  follows 
LNLN  C„  --  (mti)LN  (C„-  ^)+LM£^-  ^LM(To  (75) 

Two  axes  (x  and  y)  can  now  be  established  where 

y=  LNLN  + LKJ  CTr,  (76) 

and 

X=  IM  (G^-  Gu)  (77) 

The  correct  value  of  au  is  obtained  through  a trial  and 
error  procedure.  Various  values  are  picked  until  a 
straight  line  results  from  the  plot  of  Equations  (76)  versus 
(77)-  If  au  is  too  low,  a curve  concave  down  will  result. 
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if  too  high  the  result  will  be  a curve  concave  up  [l6]. 

The  slope  of  the  resulting  straight  line  then  becomes 
(m  + 1).  Once  m and  are  established, all  values  of  Equa- 
tion (70)  are  known  except  for  aQ.  This  value  crin  now  be 
determined  by  substituting  the  knovm  values  back  into  Equa- 
tion ( 70) . 

This  procedure  was  also  employed  by  Daniel  and  Weil 
when  their  analytical  method  proved  impractical.  One  draw- 
back found  during  their  experimentation  showed  that  gener- 
ally a very  large  number  of  guesses  are  required  to  finally 
obtain  the  desired  linear  result.  Consequently  they  de- 
vised a process  in  which  the  upper  and  lower  end  slopes  of 
the  concave  curves  resulting  from  several  estimates  of  a ^ 
are  plotted  against  a^.  This  procedure  is  illustrated  in 
Figure  9.  Connecting  the  end  slope  points  with  a straight 
line  results  in  a plot  as  shown  in  Figure  9d.  The  inter- 
secting point  provides  the  desired  values  for  m and  au« 

This  procedure  significantly  reduces  the  total  number  of 
estimates  to  identify  a^. 

The  development  of  the  methodology  for  the  parameter 
determination  for  this  study  was  based  on  two  sets  of  ex- 
perimental data,  each  with  established  material  parameters. 
One  data  set  was  for  a ceramic  material  and  the  other  for 
a composite  material. 

The  least  squares  analytical  method  of  parameter  deter- 
mination was  not  investigated  due  to  the  findings  by  Daniel 
and  Weil  [7].  A detailed  analysis,  however,  was  conducted 
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Pig.  9.  Parameter  Determination  Method 


j 
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to  establish  an  effective  way  to  use  the  graphical  approach. 
The  ceramic  material  fracture  data  was  used  exclusively  in 
this  effort  due  to  its  more  manageable  magnitudes  in  stress 
values.  The  composite  data  allowed  for  additional  verifi- 
cation of  the  final  procedures. 

The  basic  approach  of  the  analysis  followed  the  proce- 
dures outlined  by  Weibull,  with  modifications  made  as  re- 
quired. The  following  paragraphs  outline  the  sequence  of 
events  and  the  rationale  which  led  to  the  final  methodology 
for  determining  the  material  parameters. 

From  the  outset,  of  the  analysis,  a variety  of  values 
were  plotted  using  Equations  (?6)  and  (77)  as  plot  coordi- 
nates. The  resulting  graphical  data,  however,  showed  re- 
finement to  the  method  was  required.  For  all  values  of  au 
chosen,  no  singularly  definable  straight  line  or  identifiable 
concavity  could  be  established.  Figure  10  shows  a typical 
case  of  the  experimental  data  when  plotted  using  Equations 
(76)  and  (77). 

To  improve  the  method,  a smoothed  curve  was  fitted  to 
the  experimental  data.  Plotting  the  values  along  this 
curve  instead  of  the  individual  data  points,  resulted  in  a 
definably  straight  line  from  which  a ^ and  m could  be  obtain- 
ed. The  problem  with  this  technique,  however,  is  the  arbi- 
trary approach  used  in  fitting  a smooth  line  to  the  experi- 
mental data.  This  became  more  evident  when  the  published 
material  parameters  for  the  data  could  not  be  duplicated. 

A more  regimented  approach  for  data  smoothing  was  in- 
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Pig.  10.  Typical  Example  of 
Weibull  Parameter  Determination  Plot 
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stituted  through  the  use  of  the  least  squares  principle  and 
thus  various  degree  polynomial  equations  were  studied  which 
adequately  approximated  the  data.  This  required  a decision 
as  to  which  degree  polynomial  provided  the  best  results.  A 
mathematical  best  fit  could  not  be  used  since  this  would 
necessarily  be  a high  degree  polynomial,  resulting  therefore 
in  an  irregular  curve  rather  than  the  desired  smooth  curve. 
Rather  than  limit  the  flexibility  of  this  method,  numerous 
polynomials  were  incorporated  to  approximate  the  data,  and 
then  evaluated  with  respect  to  the  vaiability  in  the  Weibull 
parameters.  Although  all  higher  degree  polynomials  produced 
distinctive  values  for  m and  each  polynomial  resulted 
in  distinctively  different  values.  The  published  parameters 
were  again  not  duplicated. 

Analysis  of  the  characteristics  of  each  polynomial 
curve  revealed  an  invariance  in  the  curves  in  the  region 
beyond  the  first  three  to  six  data  points.  A very  high 
variability  was  found,  however,  in  fitting  these  first  few 
points,  figure  11  depicts  this  area.  The  characteristics 
of  smoothing  within  this  zone  resulted  in  the  following;  a 
smoothed  curve  which  flared  as  shown  by  curve  1 in  Figure 
11,  resulted  in  ou  values  near  zero  or  even  negative,  curves 
within  the  shaded  area  resulted  in  values  falling  between 
zero  and  the  lowest  failure  data  point  depending  on  whether 
the  curve  fit  was  chosen  in  the  upper  or  lower  portion  of 
this  shaded  section,  finally  if  a fit  approximated  curve  2, 
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values  of  au  greater  than  the  lowest  failure  data  point 
were  required  to  obtain  a straight  line  plot  using  Equations 
(76)  and  (77) • As  a result  of  these  characteristics  and  due 
to  its  physical  significance  within  the  Weibull  equation, 
au  was  chosen  as  a bounding  value.  The  lower  limit  of 
was  chosen  to  be  zero  based  on  the  physical  interpretation 
that  stress  below  this  value  must  result  in  a probability  of 
fracture  equal  to  zero.  The  lowest  experimental  failure 
stress  determined  the  upper  limit  for  au»  thereby  specify- 
ing that  a finite  probability  of  fracture  exists  above  this 
value. 

Using  various  degree  least  square  polynomials  a variety 
of  0u's  and  m's  were  established  within  these  limits  of  au* 
If  these  values  are  displayed  in  graphical  form  a linear 
relationship  results  as  shown  in  Figure  12.  This  data  im- 
plied that  an  infinite  number  of  combinations  of  both  m and 
0u  existed  within  the  set  limits,  any  of  which  could  be  used 
to  theoretically  duplicate  the  experimental  data.  The  pub- 
lished material  parameters  for  the  sample  data  were  also 
located  on  this  line.  Based  on  these  interesting  results, 
this  technique  had  provided  the  means  of  duplicating  the 
published  material  parameters  a3  one  of  a very  large  set. 

To  verify  that  each  set  of  values  resulted  in  a valid 
theoretical  prediction  of  the  experimental  distribution, 
randomly  picked  parameter  sets  were  plotted  using  the  Wei- 
bull equation  and  overlayed  on  the  experimental  data.  Vis- 
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ually  adequate  fit  was  achieved  with  all  sets.  However,  it 

was  noted  that  sets  with  values  of  a near  zero  or  near  the 

lowest  fracture  data  point  did  not  match  the  experimental 

data  as  well  as  those  with  intermediate  a values.  It  ap- 

u 

peared  therefore  that  a "best"  set  of  parameters  existed 
which  could  be  mathematically  isolated. 

The  mathematical  analysis  to  isolate  the  best  set  of 
parameters  employed  the  least  squares  method  of  summing  the 
squares  of  the  mean  difference  between  the  experimental  data 
and  theoretical  Weibull  equations.  The  error  between  these 
two  models  was  calculated  and  plotted  for  a variety  of 
material  parameter  sets.  The  dashed  line  shown  in  Figure 
12  represents  this  sum  of  the  mean  difference  squared  plot- 
ted against  the  corresponding  m and  a u values.  As  can  be 
seen,  a ’’best"  set  of  parameters  does  exist  among  the  family 
of  sets  represented  by  the  straight  line  when  analyzed  in 
this  manner.  Eventhough  the  set  which  resulted  from  this 
procedure  did  not  match  the  published  set  of  parameters, 
the  magnitude  of  the  overall  error  in  fit  difference  is 
very  small.  Figures  13  and  14  provide  a graphical  repre- 
sentation of  the  experimental  data  points  overlayed  using 
the  published  material  parameters  (Figure  13)  and  the  best 
set  determined  by  the  aforementioned  procedure. 

The  final  step  in  establishing  the  procedures  for  the 
parameter  determination  was  to  check  the  sensitivity  of  fit 
accuracy  between  the  experimental  data  and  theoretical  curve 
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Fig.  13 • Comparison  of  Sample  Experimental  Data  and 
Theoretical  Cumulative  Distribution  Function 
Using  Published  Parameters 
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0.97  O.  SO  0.02 


8 


FRACTURE  STRESS  »LO'(psi) 


Fig.  14.  Comparison  of  Sample  Experimental  Data  and 
Theoretical  Cumulative  Distribution  Function 
Using  Derived  Parameters 


due  to  small  changes  in  each  parameter  of  the  Weibull  dis- 
tribution function.  The  results  showed  that  all  changes  to 
au  and  m would  cause  an  increase  in  fit  error.  However,  if 

a is  made  the  variable  at  a fixed  m and  a the  author  was 
o u 

able  to  find  values  which  decreased  the  fit  error.  The 
cause  cf  this  lies  in  the  methodology  used  in  the  deriva- 
tion of  this  parameter.  The  procedure,  which  is  used  to 
find  aQ»  evaluates  Equation  (75)  for  each  stress/probability 
value  and  establishes  as  many  0q  values  as  there  are  data 
points.  These  values  of  aQ  are  then  averaged  to  obtain  the 
final  aQ-  This  final  number  therefore  is  very  data  sensi- 
tive since  the  addition  of  another  failure  value  would  gen- 
erally result  in  a change  to  aQ« 

Based  on  this  information,  the  value  of  a for  each 

o 

3e t of  m and  au  was  allowed  to  vary  approximately  plus  and 
minus  one  percent,  and  a new  fit  error  was  calculated  for 
each  resulting  new  set  of  parameters.  Although  the  reduc- 
tion in  error  is  extremely  small  in  overall  magnitude 
(.02655698)  versus  (.02655657)  the  resulting  set  of  "best" 
parameters  i3  changed  significantly  (m  = 4.868,  au  - 1661,5 
0 - 2342.06  versus  m = 4.856,  au  = 1670,  oQ  ~ 2333.52). 

Due  to  this  sensitivity  found  in  the  aQ  parameter,  the 
procedure  of  graphically  establishing  m and  au  was  revised. 
The  writing  of  a computer  program  became  necessary,  incorpo- 
rating the  experimental  data  in  conjunction  with  Equation 
(?6)  and  (77),  to  allow  a numerical  refinement  for  various 
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values  of  a . From  this  numerical  data, initial  values  of 

m could  then  be  obtained.  For  each  au* the  values  of  m and 

a are  allowed  to  vary  until  a combination  of  the  three 
o 

values  is  found  which  provides  a minimum  fit  error  to  the 
experimental  data.  Figure  15  is  a three  dimensional  repre- 
sentation of  the  physical  interpretation  of  this  error  min- 
imization process.  This  procedure  greatly  simplifies  the 
more  involved  and  subjective  graphical  method.  It  also  pro- 
vides immediate  mathematical  comparisons.  The  computeriza- 
tion verified  the  high  sensitivity  of  fit  accuracy  to  vari- 
ations in  oQ  when  establishing  an  accurate  theoretical  model 
for  an  experimental  failure  distribution. 

Since  the  above  methodology  demonstrated  that  for 
given  values  of  au,  the  fit  error  could  be  minimized,  the 
analytical  method  of  maximums  and  minimums  should  also  be 
able  to  predict  this  point.  Therefore  an  analytical  check 
of  these  procedures  was  also  conducted.  The  minimum  error 
expression  can  be  written  as 

£ £ (s-5n)z=  MINIMUM  (78) 

were 


n 


N + 


(79) 


and  Equations  (1)  and  (6). 


Taking  the  derivative  of  this  equation  with  respect  to  the 
variables  yields  the  following  expressions 
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(80) 


[(s-Sn^Sn  *0 

ni 

£(s-S")3SV*  ' 0 

n*<  ' 

l (S-S n^%-0  O 

Convergence  of  these  equations  was  verified  for  all  param- 
eters evaluated. 

The  computerized  analysis  method  which  resulted  from 
the  refinement  of  the  graphical  trial  and  error  procedure 
is  used  in  analysing  the  data  gathered  for  this  thesis. 

The  methodology  is  mathematically  reproducible  for  a fixed 
set  of  experimental  data,  and  therefore  less  subject  to 
many  inherent  inaccuracies  of  graphical  analysis  techniques. 
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IV.  Results 


This  section  contains  a compilation  of  results  from 
the  various  analyses  and  data  evaluations  comprising  this 
thesis.  The  discussion  is  divided  into  three  subsections. 

The  first  two  parts  present  the  results  obtained  from  the 
Weibull  parameter  determination  analysis  and  the  theoretical 
equation  validation  using  the  experimentally  obtained  fail- 
ure distribution.  Peculiarities,  which  are  reported  in  the 
first  two  efforts,  led  to  a supplemental  data  analysis  to 
further  investigate  the  usability  of  the  Weibull  theory. 

The  findings  of  this  evaluation  are  discussed  in  the  last 
part  of  this  section.  The  detailed  data  from  the  experi- 
mental phase  of  this  study  are  presented  in  tabular  form  in 
Appendix  B. 

Due  to  the  complexity  in  describing  the  large  assort- 
ment of  analysis  variations  which  occurred  during  the  con- 
duct of  this  thesis,  the  author  feels  it  is  necessary  to 
present  the  reader  with  an  overview  of  the  total  effort  at 
the  outset  of  this  section.  This  is  presented  in  Table  1. 
Also,  a factor  which  must  be  pointed  out  at  this  time  is  a 
peculiarity  which  was  discovered  during  the  analysis  of  the 
90-degree  flexure  specimen.  It  was  discovered  that  all 
specimen  which  had  been  equipped  with  strain  gages  had  fail- 
ed at  an  above  average  stress.  Consequently,  it  was  decided 
to  analyze  this  data  as  comprised  of  two  independent  samples* 
one  containing  21  samples,  the  other  containg  26.  All  90- 
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♦ 26  Specimen  Data  Sample 
**  21  Specimen  Data  Sample 
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degree  bending  data  will  be  analized  in  this  manner. 
Parameter  Determination 

The  flexure  data  shown  in  Tables  IX  and  X of  Appendix 
B is  used  for  the  basic  parameter  determination.  The  com- 
puterized methodology  described  in  Section  III  was  employed 
in  this  analysis.  Figures  16,  18,  and  20  present  the  final 
sets  of  m and  values  for  the  0 degree  and  two  sets  of  90- 
degree  specimen  respectively.  Additionally,  Figures  1?»  19 
and  21  provide  the  values  of  aQ  which  correspond  to  these 
values.  Superimposed  on  Figures  16,  18  and  20  are  the  fit 
error  values  which  result  from  the  correlation  between  the 
experimental  and  theoretical  distributions  when  using  the 

sets  of  m,  a and  a . Table  II  is  a summary  of  the  best 
u o 

parameter  values  obtained  as  well  as  comparison  cf  the  ex- 
perimental mean  compared  to  the  theoretical  mean  using  the 
respective  parameters. 


Table  II . 

Weibull  Parameters  Using  Bending 


Specimen 

Type 

m 

0U 

psi 

Oo 

psi 

am  ps^ 

Experimental/ 

Theoretical 

0-degree 

18.70 

0 

180010.26 

249405-50 

249307.14 

90-degree« 

.90 

8498 

1.30 

9792.25 

9736.34 

90-degree** 

.91 

8175 

2.95 

10173.09 

9981.96 

**  26  Specimen  Data  Sample 


5** 


u 


* 21  Specimen  Data  Sample 


G-degree  Pure  Bending  Parameters 
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Fig.  20.  90-degree  Pure  Bending  Parameters 
( including  strain  gaged  specimen) 


Fig.  21.  90-degree  Pure  Bending  Parameters 
( including  strain  gaged  specimen) 


These  values  are  determined  based  on  a fit  error  conver- 


gence factor  of  .0001  using  the  least  squares  method.  Table 
III . gives  the  fit  error  values  corresponding  to  the  param- 
eters shown  in  Table  II. 


Table  III. 

Bending  Data  Fit  Error 


Specimen 

Fit 

Type 

Error 

0 -degree 

.02198 

90-degree 

.04589 

90-degree 

.06114 

The  author  questioned  the  difference  in  fit  resulting  be- 
tween the  21  specimen  and  26  specimen  90-degree  parameters, 
but  one  has  to  realize  the  relationship  of  the  parameters 
in  the  Weibull  equation  and  the  significant  change  in  the 
shape  of  the  distribution  curve  caused  by  the  addition  of 
the  5 instrumented  specimen  to  the  top  of  the  distribution. 

Computations  became  some-what  erratic  as  estimates  of 
au  approached  the  lowest  failure  data  point  where  probabi- 
lity of  failure  values  predicted  by  the  Weibull  expression 
became  very  sensitive  to  very  small  variations  \r-  the  in- 
dividual parameters.  The  author  has  attempted  to  attain  as 
high  an  accuracy  as  possible,  but  practical  considerations 
required  the  establishing  of  tolerances.  As  a consequence 
of  this,  the  following  minimums  where  placed  on  the  param- 
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eter  variation  increments*  °u  t 1.0  psi j m + . 01 j qq  + .1 
psi  and  +1.0  psi.  Two  minimums  were  placed  on  a due  to 

s 

the  magnitudes  of  this  parameter  for  the  0-degree  specimen. 
Error  variation  per  unit  change  in  aQ  for  the  0-degree  data 
was  less  than  .00001. 

Figures  22,  23  and  24  show  the  experimental  pure  bend- 
ing data  and  the  Weibull  equation  plotted  using  the  values 
shown  in  Table  II.  For  comparative  purposes  additional  plots 
with  various  parameter  sets  chosen  from  Figures  16,  18,  and 
20  are  provided  in  Figures  1C  through  8C  of  Appendix  C. 

The  determining  factor  for  the  best  0-degree  specimen 
parameters  is  the  lower  limit  previously  placed  on  a . 
Mathematically,  however,  a better  fit  can  be  achieved  using 
negative  values  of  au«  The  minimum  fit  error  occured  at 
au  = -290,000  psi.  This  value  is  approximate  since  a de- 
tailed analysis  into  establishing  this  exact  value  was  not 
conducted.  Figure  2C  in  Appendix  C verifies  that  using 
negative  values  does  provide  a mathematically  accurate  fit 
to  the  experimental  data.  Within  the  physical  interpreta- 
tion of  this  parameter  as  the  zero  probability  failure 
value,  this  of  course  cannot  occur.  The  reason  for  the 
increased  accuracy  using  a negative  au  value  is  again  due 
to  the  mathematical  relationship  of  the  three  parameters 
within  the  Weibull  bending  equation.  This  of  course  also 
points  to  the  fact  that  placing  a physical  significance  onto 
the  cu  parameter  a3  it  appears  in  the  Weibull  equation,  may 
be  incorrect. 
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PROBABILITY  OF  FRACTURE 


Fig.  22.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen 


raoARaiL inr  of  frrcture 
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Fig.  23.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 

* without  strain  gaged  specimen 
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Fig.  24.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 


* including  strain  gaged  specimens 
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Theoretical  Equation  Validation 


The  parameters  developed  in  the  previous  section  now 
permit  the  evaluation  of  the  Weibull  equations  which  where 
derived  for  tension  and  three-point  loading  failure.  Based 
on  the  premise  that  the  Weibull  parameters  are  material  re- 
lated values,  the  resulting  theoretical  distributions  should 
approximate  the  distribution  established  experimentally  for 
each  failure  mode.  The  specific  evaluation  results  are  pre- 
sented per  failure  criteria  and  specimen  fiber  orientation 
in  the  following  subsect5 ons. 

A.  Tension  equation  - 0-degree  specimen. 

Equation  (15)  was  evaluated  in  combination  with  the 
material  parameters  derived  from  the  0-degree  bending  fail- 
ure distribution  as  given  in  Table  II.  Direct  substitution 
of  this  parameter  set  into  Equation  (15)  resulted  in  a fail- 
ure distribution  predicting  100  percent  failure  at  stress 
values  above  230,000  psi.  This  of  course  does  not  match 
the  experimental  data  for  0-degree  tension  failure  as  can 
be  seen  in  Table  XI  of  Appendix  B. 

Since  the  derivation  of  the  Weibull  equations  assumed 
a volumetric  flaw  distribution,  the  volumetric  differences 
between  the  bending  and  tensile  specimen  were  normalized 
using  the  following  relationship 

Vb  Vt 

a?  * a?  <81) 

where  b - bending 
t - tension 
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from  which  a volume trically  normalized  oq  can  be  defined  as 


Cob  Vt 


im 


(82) 


Equation  (15)  can  now  be  written  as 


Using  this  equation  and  the  parameter  set  established  from 

'I't 

0-degree  bending  specimen* the  resulting  theoretical  distri- 
bution function  improved  slightly,  however,  the  overall 
data  match  was  still  very  poor. 

Comparison  of  the  experimental  failure  distributions 
obtained  for  bending  and  tension  indicated  that  the  basic 
failure  stress  extremes  are  nearly  the  same.  Therefore, 
based  on  this  characteristic,  a theoretical  curve  similar 
to  that  obtained  for  bending  should  also  provide  very  ade- 
quate predictions  for  the  experimental  tension  distribution. 
The  modification  of  the  tension  equation  under  these  cir- 
cumstances can  now  be  accomplished  as  follows i since  the 
parameter  set  being  evaluated  has  a au  value  equal  to  zero, 
the  bending  risk  of  rupture  equation  can  be  written 


(84) 


The  tension  risk  of  rupture  equation  under  the  same  condi- 
tions Is 


B,-  v4  («" 


(85) 
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Similar  risk  of  rupture  values  are  obtained  for  a given 
stress  level  only  under  the  conditon  where 


(86) 


therefore,  defining 

<r.r  ■ [^5  S.  (87) 

» 

0A+  is  used  to  differentiate  this  cr  from  that  used 
o t o 

in  the  expression  normalizing  the  volumes.  By  substituting 
this  expression  into  the  tension  risk  of  rupture  equation, 
the  resulting  form  of  this  equation  is 


B. -*(•£)"  (88) 

for  the  specific  set  of  experimental  data  being  evaluated. 
This  equation  yields  very  good  failure  predictions  for  the 
experimental  data,  as  is  shown  in  Figure  25.  For  graphite- 
epoxy  specimen  with  O-dgree  fiber  orientation  it  is  there- 
fore possible  to  predict  the  probability  of  failure  in  ten- 
sion using  a set  of  Weibull  parameters  which  have  been  es- 
tablished using  experimental  flexure  data.  However,  the 
0Q  parameter  must  be  modified  to  account  for  the  differences 
in  the  nature  of  the  two  types  of  risk  of  rupture  equations. 
The  method  utilized  to  modify  the  tension  equation  has  also 
been  used  by  Hahn  and  Knight  [12]  in  studying  short  fiber 
composites. 

B.  Tension  equation  - 90-degree  specimen. 

Equation  (15)  was  evaluated  using  both  the  set  of  ma- 
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PROBABILITY  OP  FRACTURE 


i 1 — -n 1 1 i r 

>0.00  210.00  220.00  290.00  240.00  260.00  260.00  270.00 

FRACTURE  STRESS  ■ 10* ( psi ) 


280.00 


Fig.  25.  Comparison  of  Experimental  Tension  Data  and 
Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen  - using  Bending  Parameters 


terial  parameters  d<  fined  using  the  90-degree  flexure  data 
which  included  the  strain-gaged  specimen  data  points  as 
well  as  the  set  which  was  established  using  only  noninstru- 
mented  specimen  data.  Even  after  manipulations  of  Equa- 
tion (15)  in  the  manner  described  for  the  O-degree  tension 
data, no  approximation  of  realistic  failure  probabilities 
could  be  obtained  using  the  parameters  from  either  flexure 
data.  The  immediately  obvious  reasons  for  this  are  the 

values  of  a . These  values  are  greater  than  bO  percent  of 
u 

the  experimental  data  points.  As  a result, a negative  ex- 
ponentiated term  appears  in  the  risk  of  rupture  equation. 

Various  combinations  of  Woibull  parameters  with  0^ 
values  lower  than  the  lowest  failure  stress  data  point.,  ob- 
tained using  Figures  and  -0,  were  analysed  but  no  suit- 
able combination  could  be  established.  Due  to  tin?  signif- 
icant differences  in  the  failure  characteristics  in  the  00- 
degrt e fiber  oriented  graph! te- epoxy  specimen  when  subject- 
ed to  tension  as  compared  to  pure  bending  failure,  no  cor- 
relation was  found  in  the  use  of  a single  set  of  Woibull 
parameters  for  failure  probability  prediction  across  these 
failure  modes. 

C.  Throe-point  loading  equation  - 0-dogroe  specimen. 

Equation  (no)  was  evaluated  using  the  WeibuLl  param- 
eters defined  using  the  0 -degree  experimental  pure  bending 
data.  Since  the  value  of  ou  iu  this  set  of  parameters  is 
equal  to  zero,  Equation  (00)  can  be  rewritten  in  a much 
simpler  form  as  follows 
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m 


(89) 


where  a _ . is  normalized  for  volume. 
o3pt 

Substitution  of  the  parameters  shown  in  Table  II  into 
this  form  of  the  risk  of  rupture  equation  yielded  failure 
predictions  far  too  low.  For  example,  replacing  a with  the 
highest  experimentally  obtained  failure  stress  for  this 
failure  mode  resulted  in  a probability  of  failure  of  less 
than  .22.  Manipulation  of  Equation  (89)  in  the  manner  de- 
scribed for  the  evaluation  of  the  tension  equation , allows 
the  adjustment  of  the  ao  parameter  as  follows 


(90) 


(91) 


using  this  value  for  in  the  three-point  loading  risk  of 
rupture  equation  provides  the  best  correlation  attainable 
for  this  set  of  parameters.  The  resulting  curve  is  shown 
in  Figure  26.  As  is  evident, this  curve  provides  a poor  cor- 
relation to  the  experimental  data  and  therefore^  it  is  con- 
cluded that  theoretical  prediction  of  a 0-degree  three-point 
loading  failure  mode  using  0-degree  bending  parameters  is 
not  possible. 

D.  Three-point  loading  equation  - 90-degree  specimen. 

Evaluation  of  Equation  (66)  was  attempted  using  both 


I 
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— ME l BULL  CURVE 
0 S=N*/N+i 

3 - FRACTURE  PROBABILITY 
N - NUMBER  OF  3PEC1RENS 
N«  - SPECIMEN  NUMBER 

m = 18.7 

au=  0 psi 

a0=  159151.19  psi 


koo.oo  tio.oo 


tCO.OO  290.00  240.00 

FRACTURE  5TRECG 


260.00  260.00 
*10*(  psi) 


270.00  260.00 


Fig.  26.  Comparison  of  Experimental  Three-point  Loading 
Data  and  Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen  - using  Bending  Parameters 


sets  of  90-degree  bending  parameters.  Due  to  the  complexity 
of  the  basic  form  of  Equation  (66),  no  solution  could  be  ob- 
tained using  the  values  as  listed  in  Table  n.  However, 
since  the  criteria  for  establishing  Figures  18  and  20  is 

that  all  sets  of  the  listed  m,  a and  g values  provide  a 

u o 

fit  with  the  corresponding  error,  then  utilizing  the  set  of 
values  corresponding  to  qu  = 0 should  provide  an  adequate 
measure  of  whether  correlation  between  the  experimental  90- 
degree  three-point  loading  and  the  three-point  loading  equa- 
tion given  in  Equation  ( 89)  exist.  The  resulting  theoret- 
ical curves  are  shown  in  Figures  2?  and  28.  As  can  be  seen 
in  Figure  28,  the  theoretical  prediction  curve  using  the 
parameters  derived  from  the  90-degree  flexure  data  without 
the  instrumented  specimen  is  fairly  good.  However,  the  use- 
ability  of  the  theoretical  data  in  Figure  2?  is  questionable. 
The  reader  should  be  aware  that  these  curves  were  plotted 
by  direct  substitution  of  the  material  parameters  from  the 
flexure  analysis  into  Equation  (89)  using  Equation  (S?). 
Manipulation  of  Equation  (89)  as  shown  in  Equation  (9l)  was 
not  necessary.  This  of  course  means  that  the  parameters 
are  true  material  parameters  for  the  90-degree  bending  to 
90-degree  three-point  data.  The  reader  should  recall  that 
for  the  0-degree  data,  the  nature  of  the  failure  distribu- 
tion had  to  be  considered  in  order  to  achieve  correlation. 

A point  should  be  made  at  this  time,  and  that  is  that  the 
three  bending  derived  material  parameters  are  not  common  to 
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PROBABILITY  OF  FRACTURE 


**40.00  910.00  990.00  1060.00  1120.00  i 190.00  1260.00  1390.00  1400.00 

FRACTURE  STRESS  »10* (psi) 


Fig.  27,  Comparison  of  Experimental  Three-point  Loading 
Data  and  Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen  - using  Bending  Parameters* 

* 90-degree  data  ~ including  strain  gaged  specimen-. 
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PROBABILITY  OF  FRACTURE 

0.18  0.88  0.97  0.50  0.S2 


Fig.  28.  Comparison  o f Experimental  Three-point  Loading 
Data  and  Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen  - using  Bending  Parame  ters* 

* 90-degree  data  - without  strain  gaged  specimens 
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this  group  of  experimental  data,  but  are  more  oriented  to 
the  type  of  test.  As  a summary  of  this  section  the  follow- 
ing correlation  overview  is  provided. 


Table  IV. 

Correlation  Results  - Bending  Parameters 


Equation 

Data 

Fit  Obtained 

Bending 

0-degree 

Excellent 

90-degree* 

Excellent 

90-degree** 

Excellent 

Tension 

0-degree 

Good 

90-degree 

None 

Three- 

Point 

0-degree 

None 

90-degree 

Poor  using  90  flexure* 

Good  using  90  flexure**  , 

**  not  including  strain-gaged  specimens. 

# including  strain-gaged  specimens. 

Alternate  Parameter  Determination  Methods 

The  problems  encountered  in  obtaining  correlation  be- 
tween the  experimental  data  and  theoretical  predictions 
using  bending  parameters  led  to  the  evaluation  of  using  an 
alternate  parameter  determination  method.  Since  Equation 
(15)  could  easily  be  rewritten  in  the  form  of  Equation  (75)? 
it  could  therefore  be  used  in  the  same  manner  as  the  bend- 
ing equation  for  parameter  determination.  The  total  param- 
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eter  determination  effort  was  therefore  re accompli shed  us- 
ing the  tension  data.  Figures  29  and  31  provide  the  values 
for  m and  au  while  Figures  30  and  32  provide  the  values  for 
which  resulted  from  this  analysis.  Table  V provides  a 
listing  of  the  best  fit  values  contained  on  these  charts. 
Table  VI  lists  the  fit  error  corresponding  to  the  parameters 
shown  in  Table  V. 


Table  V. 

Weibull  Parameters  Using  Tension 


Specimen 

Type 

m 

a 

u 

psi 

a 

0 

psi 

0 psi 

Experimental/ 

Theoretical 

0-degree 

17.53 

0 

229965.27 

251869.99 

252128.76 

90 -degree 

4.95 

2120 

3604.64 

7228.64 

7222.43 

Table  VI. 

Tension  Data  Fit  Error 


Specimen 

Type 

Fit 

Error 

0-degree 

90-degree 

.03782 

.03316 

Figures  33  and  3^  show  the  resulting  theoretical  curves 
overlayed  on  the  0-degree  and  90-degree  specimen  tension 
data  points.  Evaluation  of  parameter  correlation  to  the 
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Fig*  33«  Comparison  of  Experimental  Tension  Data  and 
Theoretical  Cumulative  Distribution  Function 
0-degree  Specimen 
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— HE [BULL  CURVE 
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8 - FRACTURE  PROBABILITY 
H - RUflBER  OF  3PECIRENS 

mi  - spec  men  umber 

m = 4.95 
a - 2120  psi 
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Pig*  34.  Comparison  of  Experimental  Tension  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen 


other  failure  modes  for  each  specimen  type  is  presented  be- 
low. 

A.  Bending  equation  - 0-degree  specimen. 

Equation  (6)  was  evaluated  using  the  0-degree  param- 
eters presented  in  Table  V.  Results  similar  to  those  ob- 
tained for  the  previous  discussed  analysis  correlating  bend- 
ing to  tension  using  the  bending  parameters  were  obtained 
in  this  effort.  Using  the  definition  provided  in  Equation 
(87)  and  solving  the  equation  for  a the  following  express- 
ion results 

r . — /-g-T 

Db  2(rr\+')  [do*J  (92) 

The  theoretical  curve  obtained  using  this  expression  is 
shown  in  Figure  35.  This  result  show  that  parameters  de- 
rived from  either  tension  or  bending  yield  favorable  com- 
parisons to  either  theoretical  equation. 

B.  Bending  equation  - 90-degree  specimen. 

Both  of  the  90-degree  bending  data  samples  were  eval- 
uated in  this  analysis.  Poor  correlation  resulted  for  all 
equation  manipulations  attempted.  This  is  further  substan- 
tiated by  the  fact  that  the  parameter  sets  in  Figures  18, 
and  20  compared  to  those  in  Figure  31  for  a given  value  of 
<JU  do  not  agree.  The  resulting  best  theoretical  curves 
attained  are  shown  in  Figures  36  and  37. 

C.  Three-point  loading  equation  - 0-degree  specimen. 

Equation  (66)  in  the  form  of  Equation  (89)  was  evalu- 


ated.  No  correlation  was  obtained  using  the  O-degree  ten- 
sion parameters  from  Table  V.  This  is  expected  based  on 
the  results  of  the  previous  analysis  using  the  bending  pa- 
rameters. A distribution  with  a very  small  range  of  maximum 
to  minimum  failure  stress  values  would  be  expected  to  re- 
quire a high  value  for  the  flaw  density  exponent  (m)  in 
order  to  establish  correlation  between  the  theoretical  curve 
and  experimental  data. 

D.  Three-point  loading  equation  - 90-degree  specimen. 

Equation  (89)  was  evaluated  using  Weibull  parameters 
corresponding  to  0^  = 0 from  Figure  31.  Poor  correlation, 
as  shown  in  Figure  38,  resulted.  This  is  expected  since  the 
90-degree  bending  and  90-degree  tension  produced  poor  cor- 
relation, while  90-degree  bending  compared  favorably  with 
90-degree  three-point  loading. 

In  summary,  the  results  obtained  in  this  section  also 
did  not  produce  a set  of  parameters  which  allowed  fracture 
prediction  through  the  range  of  desired  failure  modes.  The 
resulting  set  of  tension  data  derived  parameters  proved  to 
be  even  less  suitable  than  those  derived  from  the  bending 
data.  An  overview  correlation  summary  is  in  Table  VII. 

Since  the  general  results  of  this  thesis  were  unable 
to  produce  a singular  set  of  Weibull  parameters  from  which 
suitable  theoretical  failure  probabilities  could  be  estab- 
lished for  the  range  of  failure  modes  investigated,  the 
author  felt  that  it  would  be  interesting  to  at  least  estab- 
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Fig.  35*  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
O-degree  Specimen  - using  Tension  Parameters 
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Fig.  36.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen*  - using  Tension  Parameters 

* including  strain  gaged  specimen* 
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Fig.  3?«  Comparison  of  Kxperlmental  Pond  in#;  Data  and 
Theoretical  Cumulative  Distribution  Kune  lion 
90-degree  Specimen*  - u;  ing  Tension  Parameters 

• without  strain  gaged  specimen 
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Fig.  38.  Comparison  of  (Experimental  Three-point  Loading 
Data  and  Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen  - using  Tension  Parameters 
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Table  VII. 


Correlation  Results  - Tension  Parameters 


Equation 

Data 

Fi  t Obtained 

Bending 

0-degree 

Good 

90-degree* 

None 

90-degree** 

None 

Tension 

0-degree 

Excellent 

90-degree 

Excellent 

Three-Point 

O-aegree 

None 

90 -degree 

None 

• without  strain. 
••  with  strain. 


lish  one  set  of  parameters  for  each  experimental  failure 
distribution  which  allowed  accurate  theoretical  duplication 
of  that  distribution.  Parameter  sets  to  provide  excellent 
theoretical  curves  for  both  bending  and  tension  have  already 
been  established,  however,  no  set  for  the  thr»'0 -point  load- 
ing data  ha3  as  yet  been  determined.  Equation  is  not 

suited  to  the  manipulation  conducted  in  Equation  (75)*  how- 
ever, the  special  case  of  Equation  (00)  presented  by  Equa- 
tion (89)  can  be  rearranged  into  a similar  form.  This  then 
allows  the  determination  of  a b03t  set  of  m and  for 
0 and  as  such  represents  best  theoretical  prediction  of  the 
three-point  loading  data  obtained  in  thi3  thesis.  The  re- 
sulting values  are  shown  in  Table  VIII. 
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Table  VIII. 


Weibuil  Parameters  Using  Three-Point  Loading 


m 

0 

0 

Error 

u 

0 

0-degree 

42.9? 

0 

211?02.2? 

.023202 

90-degree 

12.01 

0 

6150.49 

.041770 

The  corresponding  distribution  functions  are  shown  in  Fig- 
ures 39  and  40. 
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Fig.  4 0.  Comparison  of  Experimental  Three -point  Loading 
Data  and  Theoretical  Cumulative  Distribution  Function 

90 -degree  Specimen 
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V.  Conclusions  and  Recommendations 


The  objectives  of  establishing  a set  of  Weibull  param- 
eter using  experimental  data  and  then  correlating  these 
parameters  to  other  data  samples  was  accomplished.  As  a 
result  of  the  methodology  used  in  the  parameter  determina- 
tion. comparisons  of  various  values  of  m,  a and  a are 

u o 

available  for  all  data  samples  analyzed.  Listed  below  are 
some  specific  conclusions  derived  from  this  thesisi 

(1)  The  graphical  procedures  for  parameter  determina- 
tion described  by  Weibull  and  modified  by  Daniel 
and  Weil  are  not  reproducible.  The  procedures 
result  in  approximations  of  the  material  param- 
eters. These  approximations  can  easily  be  re- 
placed by  parameter  sets  with  equal  failure  pre- 
diction ability. 

(2)  The  use  of  the  above  mentioned  procedure  does 
provide  a fair  initial  guess  at  a linear  distri- 
bution of  parameter  values  from  which,  after  math 
ematical  refinement,  a best  set  of  values  can  be 
isolated  through  the  use  of  a mathematical  fit 
difference  routine  like  the  least  square  method. 

(3)  The  determination  of  the  parameter  is  not  math 
ematical ly  possible  when  a best  fit  is  desired. 
The  parameter  must  be  determined  through  a trial 
and  error  process  in  which  cu  is  fixed  and  com- 
binations of  m and  0 are  allowed  to  vary.  Accu- 


rate  fit  becomes  very  sensitive  to  the  aQ  param- 
eter. 

(4)  The  overall  analysis  was  found  to  be  very  sensi- 
tive to  data  sample  size.  Large  variations  of 
data  points  in  any  specific  section  of  a distri- 
bution was  found  to  have  significant  effect  on 
all  parameters.  Random  variations  of  a few  data 
points , however,  have  little  impact  on  m and  au 

but  tend  to  affect  a • Addition  or  deletions  of 

o 

data  points  from  a distribution  also  tended  to 
alter  the  resulting  best  fit  set  of  parameters  at 
the  sample  size  being  used  in  this  thesis.  The 
linear  representation  of  m and  ou».  however,  was 
generally  insensitive  to  this.  As  a result  of 
this  characterization,  a minimum  sample  size 
should  be  determinable  for  which  fit  error  con- 
vergence starts  to  occur. 

(5)  Except  for  tne  O-degree  tension  and  90-degree 
three-point  loading  data  the  Weibull  equations 
using  a single  set  of  material  parameters  were  not 
correlatable  through  the  three  loading  conditions 
studied  for  graphite -epoxy.  Possible  causes  for 
these  results  could  be  attributable  to  the  rela- 
tively small  data  samples  used  and  possibly  poor 
90-degree  tension  data.  However,  the  nature  of 
the  discrepancies  does  not  totally  support  this. 
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(6)  A check  of  parameter  seta  other  than  those  pro- 
viding the  beat  fit  per  au  was  not  accomplished 
due  to  time  constraints.  It  is  conceivable  that 
a parameter  3et  of  poorer  accuracy,  not  falling 
on  the  linear  m versus  presentation  exists 
which  would  allow  better  correlation  between  sev- 
eral failure  modes. 

(?)  The  physical  interpretation  generally  placed  on 
au  is  felt  not  to  exist.  This  is  demonstrated  by 
the  mathematical  improvement  which  results  in  the 
0-degree  bending  data  when  negative  values  of  j 
are  used. 

(8)  Using  the  three  parameter  form  of  the  Weibull 
equation  was  found  to  provide  better  fit  than  the 
two  parameter  (au  assumed  zero)  form.  For  the 
data  used  in  this  thesis,  an  assumption  of  = 0 
would  have  resulted  in  a significantly  increased 
fit  error. 

(9)  There  are  a great  number  of  limitations  In  the 
use  of  the  Weibull  equations.  For  the  composite 
studied  and  sample  size  uaed,  the  Weibull  theory 
is  only  useful  if  the  parameters  are  derived  based 
on  specific  experimental  data.  Correlation  wa3 
established  between  O-degree  tension  and  bending 
only  after  manipulation  of  the  Weibull  equation 
based  on  the  characteristics  of  the  available  ex- 
perimental data.  Thus,  the  use  of  the  Weibull 
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equations  does  not  lead  to  a reduction  in  test 
requirements  and  the  associated  cost  savings  for 
composite  materials.  From  the  available  data 
used  in  this  study,  it  is  appropriate  to  say  that 
the  only  true  correlation  may  have  between  the 
90-degree  three-point  loading  and  90-degree  bend- 
ing without  instrumented  specimen  data  samples. 
However,  a more  conclusive  evidence  could  be  de- 
veloped with  additional  experimental  information. 

(10)  It  is  recommended  that  the  experimental  data  es- 
tablished in  this  study  be  extended  with  addit- 
ional test  data  to  investigate  the  variability 
of  the  Weibull  parameters  as  additional  data  is 
added.  Care  must  however  be  taken  to  insure 
specimen  homogeneity  in  both  geometry  and  material. 
Once  a substantial  sample  size  has  been  obtained, 
the  error  to  parameter  sensitivity  can  be  investi- 
gated by  random  deletions  or  additions  of  data 
points  within  this  sample. 


9? 
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APPENDIX  A 


Three-Point  Loading  Equation  Derivation 

This  appendix  provides  the  detailed  steps  required  to 
obtain  the  form  of  the  risk  of  rupture  equation  shown  for 
Equation  (11). 
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Starting  with  Equation  (8) 

b % x.  „ 


O 2, 


where 


v - In  L 
Xu  * C^TZ 


-1  3 C~u  h_ 

Z*  TbT 


Integrating  with  respect  to  Y yields 

B-WVW-crdxdz 


*m  2u 


Integrating  with  respect  to  X yields 


taking  out  the  constants 
r\  /llim  . 


B-? Wmwl  i(Sxz-c)“*'|;d 

In  1 XH  W 

substituting  the  integration  limits  yields 

B'zrAcc;  / 

Cm 

Defining 


I-  L*-4(,P^)~d* 


The  integration  of  this  integral  is  now  addressed  separa- 
tely in  the  following  manner. 


Defining 


<T  * 


(9A) 


then 


Is  / ~?i<rz  ~ m4'dz 


(10A) 


Integration  of  this  term  results  in  the  following 

T = L^Z-g:)^1  |k/t  <-c)  (rn-n)  A / . 


( 11A) 


Further  integration  of  the  resulting  integral  leads  to  a 
form  which  can  be  expressed  in  series  notation 


I-\l  ♦ t-rf-j  i(m-  C)«dr 


using  the  following  definitions 

m * m * oC 
N*  m 

leads  to  the  form  defined  in  Equation  (11) 


(12A) 


(13A) 


(14A) 


!•{[  W" ej.  {‘tW-cJ-dz 


(15A) 
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Substitution  of  this  final  expression  for  I into  the  risk 
of  rupture  equation  will  result  in  the  form  of  this  equa- 
tion shown  in  Equation  (11) ■ 


APPENDIX  B 


Experimental  Data 

This  appendix  provides  a compilation  of  the  experimen- 
tal data  gathered  for  this  study. 
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Table  IX . 

Four  Point  Loading  Data  - O-degree  Specimen 


Stress** 
(lb/in  ) 
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Table  X . 

Pour  Point  Loading  Data  - 90-degree  Specimen 
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Table  XI. 

Tension  Loading  Data  - Q-degree  Specimen 
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Table  XII. 

Tension  Loading  Data  - 90-degree  Specimen 
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Strain-gaged  Specimen  # Fractured  at  tabs 

Calculated  - using  a = P/A 


Table  XIII . 

Three  Point  Loading  Data  - O-degree  Specimen 
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APPENDIX  C 


Data  Plots 

This  appendix  provides  graphical  data  verifying  the 
linear  relationship  between  the  values  of  m and  a^.  Also, 
the  validity  of  using  the  sum  of  the  mean  error  technique 
for  isolation  of  the  best  set  of  values  can  easily  be  seen 
by  plot  comparisons  of  similar  data. 
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Fig.  1C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
O-degree  Specimen 


PROBABILITY  OF  FRACTURE 

0.97  0.60  0.02 
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Pig.  2C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
O-degree  Specimen 
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PROBABILITY  OF  FRACTURE 

O.l  2 0.2S  0.97  O.SO  0.02 


HElBULt  CURVE 
SrN«/N* l 
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- NUMBER  OF  SPEC IRENS 
i - SPECIMEN  NUMBER  * 


= 1.12 

u=  8365  psi 

o=  5*83  psi 
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FRACTURE  5TRE55  *10' (psi) 


Fig.  3C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 

* without  strain  gaged  specimen;. 
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Fig.  4C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 

* without  strain  gaged  specimens 
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PROBABILITY  OF  FRACTURE 


- ME [BULL  CURYE 
S=N*/N*i 


3 - FRACTURE  PROBABILITY 
N - NUMBER  OF  SPECIMENS 
N*  - SPECIMEN  NUMBER  • 


m = 1.46 
au=  8175  psi 

a0=  22,?6  psi 
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FRACTURE  STRESS  •lO'(psi) 


Fig*  5C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 


* without  strain  gaged  specimen 


PROBABILITY  Of  FRACTURE 
o.it  o.ts  o.s?  o.so  o.t* 


Pig.  6C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 

• including  strain  gaged  specimen*. 
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PROBABILITY  Of  FRACTURE 
O.l*  O.tS  0.97  0 .BO  0.0* 


Pig*  7C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 

* without  strain  gaged  specimen :> 


118 


PROBABILITY  OF  FRACTURE 


— MEtBULL  CURVE 


S=N«/N*l 

3 - FRACTURE  PROBABILITY 
N - NUMBER  OF  SPECIRENS 
Na  - SPECIMEN  NUMBER  • 


m = 10.61 

au=  2000  psi 

a0=  4375.16  p3i 


S 90.00  90S*  00  040.00  OH.  00  1010.00  t'tOO.OO  1100.00  ItlO.OO  IfYO.OO 

FRACTURE  STRESS  ■lO\Psi) 


Pig.  8C.  Comparison  of  Experimental  Bending  Data  and 
Theoretical  Cumulative  Distribution  Function 
90-degree  Specimen* 

* without  strain  gaged  specimen.: 
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